Decoherence in a superconducting flux qubit with a 7r-junction 
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We consider the use of a 7r-junction for flux qubits to realize degenerate quantum levels without 
external magnetic field. On the basis of the Caldeira-Leggett model, we derive an effective spin- 
Boson model, and study decoherece of this type of qubits. We estimate the dephasing time by using 
parameters from recent experiments of SIFS junctions, and show that high critical current and large 
subgap resistance are required for the 7r-junction to realize a long coherent time. 
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It is now well established that in addition to conven- 
tional Josephson junctions having an energy minimum 
at zero phase difference across the junction, there exist 
the so-called 7r-junctions which provide the phase shift 
of TT in the ground state. The intrinsic 7r-shifts were 
first realized in grain boundary Josephson junctions in 
d- wave superconductors [JQ- Subsequently, 7r-junctions 
have been realized in hybrid structures between high-Tc 
and low-Tc superconductors 0, [1] by injection of 
quasi-particles [4j. Recent development in fabrication of 
superconductor-ferromagnet-superconductor (SFS) junc- 
tion made it possible to obtain a 7r-junction with high 
critical current density [s]. An advantage of SFS junc- 
tions is the possibility to combine them with usual low-Tc 
superconductive circuits using conventional fabrication 
technique. 

The use of 7r-j unctions provides several new applica- 
tions. For example, the application of 7r-junctions as 
complementary devices in SFQ logic was recentljj_j)ro- 
posed ^ and realized in high-Tc-low-Tc junctions 7^. It 
is interesting that before this 'classical' application of the 
TT-junction, the use of 7r-junctions for realization of quan- 
tum two-state systems was considered 0. In this qubit 
system, the 7r-junction was used as a tt phase shifter along 
the loop instead of current biasing or external magnetic 
flux. After this proposal, remarkable progress in fabrica- 
tion, coherent control of one qubit, controllable coupling 
between qubits, and readout with high fidelity has been 
achieved in superconducting qubits [1, [13, El, [13, [l3l ■ 
Nevertheless, up to now, the use of 7r-junctions to qubits 
has not been studied experimentally. One of difficulties 
for realization may lie on the original proposal in which a 
qubit consists of complicated circuits with many Joseph- 
son junctions @. Another serious difficulty comes from 
dissipation due to quasi-particle excitation, which is un- 
avoidable in many realizations of 7r-junctions. Generally, 
qubits suffer strong decoherence by excitation in the en- 
vironment. 
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FIG. 1: Flux qubit circuits with a 7r-junction. Shunt resis- 
tance at the TT-junction is introduced for estimate of damping 
effects. 



In this paper, we consider the use of 7r-j unction for 
phase bias of flux qubits. The circuit we study is shown 
in Fig. [T] In this circuit, we need no external flux to re- 
alize degenerate quantum levels, because the phase drop 
across the three Josephson junction is adjusted as tt by 
the TT-junction with a large Josephson energy. This type 
of phase bias can avoid dephasing due to noise in exter- 
nal flux, and is frequently called as a 'quiet qubit'. In 
actual experiments, however, damping at the Tr-junction 
may cause severe decoherence on the qubit. The purpose 
of this paper is to derive the effective spin-Boson model 
describing the flux qubit with a damped 7r-junction, and 
to estimate the dephasing time by using realistic exper- 
imental parameters. We clarify the condition for long 
coherence time in this qubit system, and discuss the pos- 
sibility of the use of Tr-junctions for qubits by referring 
recent experiments on SIFS junctions. 

In order to describe damped dynamics, we introduce 
the RSJ model for the 7r-junction as shown in Fig. [2 
where dissipation is expressed by a resistance R shunted 
in parallel to the junction. We expect that this phe- 
nomenological model may give a qualitative estimate of 
decoherence effects by Tr-junctions. We introduce the 
charging energy Ec,-k = e^/(2C7r) and damping fre- 
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quency 7 = I/IRCt^) of the 7r-junction. 

The Hamiltonian consists of three parts as H — 
-ffqubit + Hl + i?7r ■ The first part iJqubit describes a flux 
qubit, and is given as p3] 



H, 



qubit 



= — i?j(c0S 01 + COS </i2 + Ck cos (/)3) 



3^ 



(1) 



where Ej is a Josephson energy, and Ec is a charging en- 
ergy. Here, (pi and rii are a phase difference and induced 
charge at the z-th junction, respectively. The area of one 
junction is reduced by the factor a, which is typically 
taken as 0.8 The second part of the Hamiltonian 

describes the inductance energy of the loop, and is given 
as 
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where = h/ (2e). Here, cf>Tr is a phase of the tt junction, 
and (pext = 27r($oxt/$o) is a phase induced by the exter- 
nal flux through the loop. By assuming small inductance 
L, the inductive part of the Hamiltonian can be treated 
as a constraint condition 

4>1 + 4>2 + 4>3 + 4>7T = 4>cyit- (3) 

The third term describes the damped 7r-junctions, and is 
expressed by the Caldeira-Leggett Hamiltonian 




4^Ec,^ni 



Cq, 



■(4) 



The damping property is determined by the spectral 
function 



(5) 



In the RSJ model, the spectral function is given as 

J(w)/A/ = 7tje-"/"% (6) 

where M = l/(8i5c,7r) is a mass of the 7r-junction, and 
LOc is a high-frequency cutoff. 

In this paper, we focus on the 'passive' use of the tt- 
junction. For this use, the Josephson energy of the tt- 
junction should be taken as sufficiently large. Hence, we 
assume Ej^t, 3> Ec,E,j, and approximate the Josephson 
energy of the 7r-junction as Ej^T^{(t)T^ — 7r)'^/2. Within this 
approximation, the phase of the 7r-junction is kept almost 
TT. The remaining dynamics around the potential minima 
is described by a damped oscillator with a eigenfrequency 
ujQ = {SEj^T^EcTiY^^/h. In the following discussion, we 
set h — 1. 

Under the condition Ej ^ Ec, which is taken for 
usual flux qubits, we can truncate the Hamiltonian -ffqubit 
into the two-level Hamiltonian as -ffqubit = ^^two-stato + 



^coupling- The flrst part iJtwo-stato = {^/2)ax + {s/2)az 
describes the qubit system, where A is a tunnehng split- 
ting, and e is a bias proportional to the external flux $ext. 
The second part, which describes the coupling between 
the qubit and the 7r-junction, is given as 



coupling 



(7) 



where Ej^.e = (1 - l/{Aa^))^/^Ej, and A4>^ = 0^ - ^. 

To simplify the Hamiltonian i/^, we change the vari- 
ables as a; = Af^/^At/)^, p — M^^/^n^. We further re- 
place the sum in the Hamiltonian of the harmonic os- 
cillators by the integral. This can be performed by re- 
placing the variables as — m}^'^Xa/{ALijy/^ , — 
m-i/2p^/(Aw)i/2, and = {Mm^y^/^Ca,/ {ALof/^ , 
where Aw is a length of one slice in the w-direction. In 
the limit Aw — > 0, we obtain 



1 , 

duj 







^uj^X^ - C^X^x + ^x^ 



.(8) 



The coefficient C^j can be related to the spectral function 
as 



J{uj)/M 



2 LO 



(9) 



The Hamiltonian of the 7r-j unction describing a 
damped oscillator can be diagonalized exactly [H, [l^ . 
In order to express the eigenmodes with the energy w, 
we introduce a canonical transformation for the opera- 
tors as 

/•OO 

X^^a{uj)x+ duj'b^,{uj)X^,, (10) 
Jo 

where the coefficients, a{uj) and bi^i(uj) are assumed to be 
real. The coefficients are chosen to satisfy the eigenmode 
equations 

/"OO 

WQa(w) + / cLu'C^'b^'{uj) — Lu'^a^Lu), (11) 



(12) 



C^'a{uj) + Lu''^b^' (cj) = uj'^b^' (uj), 



where ujq^ujq+J duj'C^i2 /w'^. Then, the Hamiltonian 
of the TT-junction can be diagonalized as 



p2 ^ 



(13) 



In order to solve the eigenmode equations, eqs. (|lip and 
(fT^ . we may follow the calculation in Fano's paper p^ . 
We only give the result for a(w) as 

= t 2ni/A 2-,^?\ -2 p( ^^2 (14) 
(7r^C4/4w^) + [uj^ - t^o - Fy^)) 



F{uj) = P duj 
Jo 



(w2 - w'2) 



(15) 
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The part of the energy renormaUzation is modified as 



F{uj) = 



dio' - 



(16) 



Here, the second term in r.h.s. can be neglected, because 
it can be shown to be 0(A/wc)- 

Thus, we obtain the new expression for the vr-junction 
as (fT3| . while the coupling term ([7]) is rewritten by the 
relation 



dLu'a{uj')X^r . 

As a result, we obtain the total Hamiltonian as 
A e 



(17) 



H 



2 2 



duja{uj)Xi^ 



(18) 



In this modified spin-Boson model, the effective spectral 
function is given by 



{^Cl/2u) 



(^C2/2c.)2 + (c^2_^2)2- 



(19) 



For the RSJ model, by substituting J{uj)/M = 
7rC^/(2ti)) = ju!, the effective spectral function is ob- 
tained for w <C Wc as 



j'^Uj'^ + (W^ - Wg)^ 



(20) 



Note that the form factor of the damped oscillator (the 
factor in the bracket) appears in the effective spectral 
function. 

By using the effective spectral function J^siuj), we es- 
timate the dephasing time of the qubit at the optimal 
point (e = 0), where long coherence time is realized by 
suppressing a linear coupling to the heat-bath. The de- 
phasing time is evaluated within the spin-Boson model 
in the form [H, |T3| 



relax 



T* 



The relaxation rate Tj-ciax is calculated as 

A 



^relax = 2Jeff(A) COth 



2fcBT 



(21) 



(22) 



On the other hand, I/T2, which is a pure dephasing 
rate due to a quadratic coupling to the heat-bath at the 
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FIG. 2: The estimated relaxation time Treiax and pure dephas- 
ing time T2 ■ The result for a 10 /im x 10 (1 fj,m x 1 fim) 
TT-junction is shown by the solid (dashed) line. 



optimal point, is calculated from the coupling strength 
r = lim^^o J{oj)/uj as [3 
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(23) 
(24) 



We estimate the dephasing time in the present flux 
qubit by using the parameters in Ref. [l^. In the ex- 
periment, the parameters are chosen as Ej/ks = 12 K, 
Ec/kB = 350 mK, A/fce = 160 mK, T = 25 mK, and 
a = 0.8. There are several candidates of 7r-junctions for 
phase bias. We have estimated dephasing time for several 
TT-junction systems, and found that only underdamped tt- 
junctions may give a sufficiently long dephasing time. 

Here, we discuss underdamped SIFS junctions by using 
the parameters in Ref. [l^ We choose a capacitance and 
a subgap resistance for unit area as c = 0.08 F/m^ and 
r„ = 3.0 X 10-"^ 17m2 respectively In Ref. [H, the 
measured critical current density is jc — 5.0 x 10^ A/m^. 
Here, we take the critical current density as a parame- 
ter, and discuss its dependence keeping constant. In 
Fig. [21 we show the relaxation time Treiax and the pure 
depasing time Tj* as a function of the critical current 
density for a 1 /im x 1 fj,m and 10 fim x 10 //m junction. 
In this estimate, the relaxation process is always domi- 
nant (Treiax ^ T2), and therefore the dephasing time 
is determined by Trqi^. We find that at the critical cur- 
rent density of Ref. 19 the dephasing time is very short, 
while long coherence time is obtained for junction with 
larger area and higher critical current density. The re- 
laxation time has a resonant structure at a low critical 
current jc = j* , where the resonant condition wq = A is 
satisfied. For jc ^ jc , the relaxation time and the pure 
dephasing time depend on the junction area A and crit- 
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ical current density jc as Treiax oc j'^A and Tj* °^ Jc^^j 
respectively. As seen in Fig. [31 if we use underdamped tt- 
junctions with large critical current (jc ~ 10^ A/m^) and 
large junction area (^4 ~ 10 /imx 10 /im), coherence time 
becomes of order of 1 ms, which is sufficiently long com- 
paring to the decoherence time limited by other sources. 
We note that when the relaxation process is dominant, 
the dephasing time is proportional to Ej^R; for realiza- 
tion of long coherence time we need to increase both the 
Josephson energy Ej^^^ and the subgap resistance R of 
the TT-junction. 

Thus, for long-time coherent operations, one has to 
improve the quality factor by changing experimental pa- 
rameters of TT-junctions. Especially important parame- 
ter is the critical current density in the 7r-state of the 
junction. In the SIFS junctions described in Ref. [l^ the 
critical current density jc ~ 5 x lO'^A/m^ in the 7r-state 
was still rather low, three orders of magnitude less than 
jc ~ 4 X lO^A/m^ in a SIS junction having the same tun- 
nel barrier. Possible reason for the strong suppression of 
the critical current is the use of diluted alloy Nij-Cui-^; 
which has rather strong disorder leading to fast decay of 
the supercurrent with increasing F-layer thickness. Since 
O-TT transition occurs at certain critical thickness of the 
F-layer, the supercurrent in the 7r-state is much smaller 
than in the 0-state. However, smallness of jc is not an 
intrinsic property of SIFS junctions. In a clean homoge- 
neous fcrromagnct the decay length may become much 
longer than the O-tt transition thickness. Recent exper- 
iments [m using NiaAl have demonstrated multiple O-tt 
transitions with only modest decay of |jc| as a function 
of the thickness of NiaAl. Therefore, choosing different 
materials for a ferromagnet layer may finally lead to in- 
creasing jc and thus to an increase of the dephasing time 
of qubits with an SIFS junctions. 

Finally, we discuss the advantage of the present phase 
bias. In usual flux qubits, external magnetic flux is 
needed to produce phase bias along the loop. In many 
experiments an external coil with a large current and 
weak coupling to the qubits has been used. However, 



this prevents one from using a superconducting shield 
which provides good shielding of qubits from external 
flux noise. On the other hand, if one uses a local bias- 
ing with a control line, noise in the current source de- 
grades the coherence of the qubit. From eq. ((22)) . the 
relaxation time due to this noise can be evaluated as 
^r^L = (2A/n2)(MJp)VZ, where Ip = 2nEj^,e/<^>o is 
a circulating current, Af is a mutual inductance, and Z 
is an impedance of the current source For ob- 

taining the dephasing time longer than 1 ms, the maxi- 
mum value of the mutual inductance M is estimated to 
be 0.03 pH for Z = 50 fl. Then, the external current 
needed for the phase bias becomes 40 mA, which is un- 
realistically high. Therefore, the use of 7r-j unctions may 
be an attractive option for individual phase biasing on 
qubits. 

Recently, another phase biasing scheme with a trapped 
flux in a superconducting loop has been proposed and 
demonstrated [1^. A possible advantage of our scheme 
using a 7r-junction is that we do not need to apply the 
large external fleld corresponding to a half flux in the 
loop ever, either globally or locally. This makes imple- 
mentation of a superconducting shield simpler. 

In summary, we proposed a simple phase bias by tt- 
junctions for flux qubits, and studied dissipation effects 
at the TT-junction. In the framework of the Caldeira- 
Leggett model, we derived the effective spectral func- 
tion of the spin-Boson model, and used it for estimate 
of the dephasing time of the proposed qubit. We showed 
that for long coherent operation both the subgap resis- 
tance and critical current of the TT-junction have to be 
increased. We expect that further improvement in qual- 
ity of TF-junctions enables us to use it for a tt phase shifter 
for flux qubits. 
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